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INTERPLANETARY MIDCOURSE GUIDANCE USING RADAR TRACKING 


AND ON-BOARD OBSERVATION DATA 

By Luigi S. Cicolani 
Ames Research Center 


SUMMARY 


This paper examines the performance of midcourse guidance and navigation 
systems for manned interplanetary missions when radar range and range -rate 
measurements and on-board theodolite observations are available as sources of 
navigation data. Results are given for two Mars missions^ one a direct flight 
to Mars and back and the other a round trip with return via Venus . 

Reasonable observation schedules to refine the miss estimate to the 
accuracy required for entry or close approach are found to consist primarily 
of radar observations for estimating the greatest part of the miss^ plus fre- 
quent on-board observations of the target planet during the approach to Mars 
or Venus ^ or radar observations during the approach to the Earth. Additional 
on-board observations during the departure from distant planets are usually 
helpful in reducing fuel requirements. 

The system performance is such that the miss is generally estimated better 
than it can be controlled. Consequently^ the work of Lawden can be applied to 
obtain an analytical solution for the a priori optimum correction schedule. 

The minimiHn fuel requirements^ the correction size^ and the dependence of 
these parameters on the number of midcourse corrections are also determined by 
the theory. Results of optimimi four -correction schedules used in the simula- 
tions on all legs of the two missions agree fairly well with the theory. 

The velocity -correction requirements for the various legs range from l8 
to 2k m/sec, and the rms miss in pericenter radius at arrival is less than 
2.8 km at Earth and Mars and 6 km at Venus. The downrange miss varies from 
2 to 22 km. 


INTRODUCTION 


In manned interplanetary missions the advantages of a self-contained 
on-board midcourse guidance and navigation system independent of communica- 
tions with Earth are evident. The performance of such systems is discussed 
in reference 1^ where it is found that navigation for the entire mission can 
be carried out reasonably well with only the on-board system. 

Nevertheless^ communications \r±ll normally be operational and both 
ground-based radar tracking and on-board observations will be available as 
navigation data sources. Ground tracking is usually superior to on-board 
observations in estimating the trajectory^ but its ability to establish the 



relative motion to the accuracy desired for manned missions during the 
approach to distant planets is limited. Consequently; the present study was 
initiated to evaluate the performance possible from the use of both tracking 
and on-board observational data. 

Midcourse guidance and navigation refers to the process of estimating and 
controlling the terminal miss at the target planet. This is carried out by a 
sequence of observations and trajectory corrections. To determine system per- 
formance it is necessary^ firsts to obtain appropriate observation and cor- 
rection schedules. The observation schedules are determined empirically to 
utilize the more effective observation type or combination of types in each 
region of the mission; sufficient observations are made to determine system 
performance near the limit of its capability. The correction schedule should 
be selected to minimize the fuel requirements and^ in this study^ a good 
approximation to the optimum correction schedule is obtained analytically 
using the method of Lawden (ref. 2). 

The discussion of system performance considers the accuracy of the miss 
estimation^ the terminal miss^ and the velocity correction requirements. 

Results are given from simulation of the system for two reference missions; both 
are round-trip missions to Mars^ one a direct trip and the other includes a 
swing past Venus during the return to Earth. The trajectory estimation is 
based on the Kalman filter theory and the results are given in terms of the 
standard deviations of the random variables of interest. 
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covariance matrix^ current uncertainty in state estimate 

covariance matrix^ current state deviation 

covariance matrix^ measurement error 
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covariance matrix, velocity-correction error 
sum of irnis midcourse velocity corrections 
time 

reference arrival time 

time of jth midcourse velocity correction 
unit vectors, defined where used 
covariance matrix, velocity correction 
velocity correction 

current state deviation from reference state 
measured value of observable 

variance, aiming error in velocity correction 

small deviation from reference value of ( ) 
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value of ( ) at injection 
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e:xpected value of ( ) 
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ANALYSIS 


The performance of the midcourse guidance and navigation system is 
evaluated by simulating its operation. The detailed analysis of this simula- 
tion has been discussed extensively in reference 1 and the reports it cites; 
thus^only a brief review of the pertinent background is necessary here. 


Simulation of the Interplanetary Guidance and Navigation Problem 

The following process is to be simulated. A vehicle is injected onto a 
reference interplanetary orbit ^ that is^ one with desired end conditions. 
Because of random injection errors the reference orbit is not obtained and 
the guidance and navigation system must estimate and control the end condi- 
tions. The estimation is carried out by processing observation data. In this 
study the data are obtained from on-board theodolite observations of various 
planets and Earth -based radar range and range -rate measurements. These obser- 
vations are subject to random measurement errors and their processing is based 
on the Kalman filter theory. Several times during the flight impulsive veloc- 
ity corrections are made to control the end conditions. These corrections are 
determined by means of a guidance law from the estimated end conditions^ and 
are subject to random execution errors. 

In the mathematical analysis of this process two assumptions are made: 

(l) all input random variables (injection^ correction^ and measurement errors) 
are gauss ian distributed with zero means ^ and (2) the orbits of interest are 
sufficiently near the reference orbit to allow linearization of the equations 
of motion. Consequently^ the output random variables (velocity corrections^ 
terminal miss^ miss estimate uncertainty^ vehicle state ^ etc.) are linear com- 
binations of the input random variables and^ therefore^ are gaussian distrib- 
uted with zero means. Thus^ the simulation need consider only the covariance 
matrices of the various quantities of interest. 

Reference traje ctories.- Two round-trip manned Mars mission trajectories^ 
including a stopover on Mars^ are used as reference trajectories. Both are 
taken from reference 1 where they are referred to as the high-speed mission 
and the Venus swing -by mission. A brief list of orbital parameters is given 
in table I. Reference 1 also gives projections of these trajectories and of 
the planetary orbits on the ecliptic plane^ which are useful in visualizing 
the relative position of the vehicle and planets at various times in the 
mission. 
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Injection and velocity-correction errors * ~ The injection errors of 
reference 1 are assimed here and are listed in table II. The errors are 
assumed to be identical at both Earth and Mars injections and also are identi- 
cal to the initial uncertainty in the estimate of the state. The miss and 
miss uncertainty at the start of the Venus -Earth leg of the Venus swing -by 
mission are also listed in table II. These are^ of course^ obtained from a 
particular simulation of the preceding Mars -Venus leg. 

For velocity-correction errors^ the model of reference 1^ appendix has 
been adopted. Briefly^ it includes proportional^ aiming^ and cutoff errors^ 
and has the following covariance matrix for the correction error: 


S = K^c + ^ (c^I - C) + = 

^ c^ 

where C is the covariance matrix of the commanded velocity correction and 
c^ its trace. The values of ^ and (variances of the proportional^ 

aiming^ and cutoff errors) are listed in table II. 

Each midcourse velocity correction attempts to null the miss as computed 
from the estimated state. Fixed time -of -arrival (FTA) guidance (see ref. 1 for 
the corresponding expression for C) will be used throughout. Miss for FTA. 
guidance is defined as the position deviation at the reference arrival time^ 
which corresponds to pericenter on the reference orbit. 


Filtering of Measurements 

All navigation information is assumed to be processed by a Kalman filter. 
Let P(t) be the covariance matrix of the uncertainty in the estimated current 
deviation state. Suppose the measurement is 

y ^ Hx + e 

where y is the measured value of the observable^ Hx is the true value of 
the observable^ and e is the random error in the measurement. The H 
matrix consists of partial derivatives of the observables with respect to the 
state variables. After an observation^ the new uncertainty covariance matrix 
is 


P' = (I - KH)P 

where K is the filter weighting matrix 

K = P#(Q + HPH'T)"^ 


and 


Q =. (ee^) 
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Qn-'board observ ations. - A nimber of on-board optical observations can be 
considered; for example; planet diameter; sextant; and theodolite measurements. 
Reference 1 indicates that planet diameter measurements are inferior and that; 
while the sextant is more convenient to operate than the theodolite; a suffi- 
cient number of sextant observations with suitably chosen stars will give 
results equivalent to those obtained from the theodolite. For these reasons 
and for computational simplicity; the on-board measurements of this study are 
restricted to theodolite observations. 


The theodolite measures the celestial latitude and longitude of the 
observed body and thereby determines the vehicle’s position deviation normal 
to the line of sight to the planet. Its H matrix and error model are dis- 
cussed in reference 1. The resulting covariance matrix of measurement errors 
iS; from reference 1; the diagonal 2x2 matrix: 


Q = 


inst 


+ b^ 



I 


where 0 is the half subtended angle and R is the distance to the observed 
body. The quantities b; and 5Rg are the random instrument error; 

planet radius uncertainty; and planet position uncertainty; respectively. 
Their rms values are listed in table II. 


Radar observations . - A network of three stations (Goldstone; Johannesburg; 
and Woomera) was assumed with the viewing station measuring range and range 
rate simultaneously. The (2 x 6) H matrix for these measurements is 


H 


T 



o 


I - UpUj. 


Vv/^ 


Ur 


where u^ is a unit vector in the direction of the line of sight from the 
station to the vehicle; ^v/s velocity of the vehicle relative to the 

viewing station; and r is the range. All quantities are evaluated on the 
reference orbit. The two columns of are the derivatives of range and 

range rate; respectively; with respect to the state variables. 


Radar measurements are affected by a variety of errors from various 
sources including white noise; colored noise (noise with finite correlation 
time); bias errors from the radar equipment; and errors from other sources; 
such aS; station location; station clock; and astronomical unit uncertainties. 
A comprehensive list of such errors in a study of guidance and navigation in 
the lunar mission is considered in reference 3* '^or simplicity; the present 
study considers only white -noise errors. The variances of these errors for 
the range and range -rate measurements are given below. These were obtained 
from reference 3 by combining the white- and colored-noise variances listed 
there. This is approximately valid because the time constant of the colored 
noise is short compared to time between observations. 
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Range: + ai^r^ (km^) 

Range rate: + biCp^r^ + bscp^r + b 3 cp^ (m/sec 

where 

cp = 1 + 0 . 03 r bi = 2x10“^^ 

ao =. 0.00525 b 2 = 6. 5x10 

ai = 0.767x10"^ b3 = 2x10"^ 

bo = 0.1 

In these expressions the dimensions of r and r are km and km/sec ^ respec- 
tively. 

The rms range measiirement error increases almost linearly with range 
outside the Earth’s sphere of influence (fig. 1). 

The range- rate measurement error depends on both the range and range 
rate. As an example^ the time history of the rms range -rate error for the 
outbound leg of the high-speed mission is plotted in figure 2. 

Velocity -correction measurement . - Accurate accelerometers are assumed to 
be on board the vehicle for measuring the executed velocity corrections. The 
accelerometer measurement errors are assumed equally likely in all directions 
with an rms value of 1 cm/sec. 

Velocity -Correct ion Schedule 

With the use of radar the question of an optimum velocity -correction 
schedule can be settled analytically from the work of Lawden (ref. 2). The 
application of his work to the present problem is discussed in detail in 
appendix A^ and a brief description of the solution is given below. 

The analytical solution is based on three assumptions : 

1. At the time of each correction the miss has been estimated better 
than it can be corrected. 

2. The terminal miss due to current velocity deviations is propor- 
tional to the time to go. 

3 . The miss generated by random errors in each correction is princi- 
pally caused by the cutoff error. 

The first assumption leads to a velocity-correction schedule independent of 
the observation schedule; consequently^ the first correction cancels almost 
all the miss due to injection errors^ and subsequent corrections account prin- 
cipally for the miss generated by random errors in the preceding correction. 
The problem then becomes mathematically tractable as a result of the second 
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and third assumptions. In the solution^ the optimum timing of each correction 
represents a con^romise between the fuel cost for subsequent corrections^ which 
decreases with a delay of the current correction because the post-correction 
miss is reduced^ and the fuel cost of the current correction^ which increases 
with a delay. 

In the current simulations the three assumptions are generally satisfied. 
In all the trajectories examined here^ radar observations are sufficiently 
effective in estimating the orbit (compared to <the errors in correcting it) 
that the first condition is generally satisfied. The second assumption 
requires that the matrix A 2 (tp^t) (a submatrix from the state transition 
matrix which relates position deviations at t^ to velocity deviations at t) 
be proportional to the time to go^ typical interplanetary orbits^ 

this assumption is satisfied for t^^-t corresponding to somewhat less than 
half an orbit. This assumption is commonly violated in the early phase of an 
interplanetary orbit. However^ no difficulty arises because the optimum tim- 
ing of the first correction is determined empirically and^ therefore^ none of 
the three assumptions need be met for the first correction. The remaining 
corrections will usually occur much later in the flighty more than halfway to 
the target position^ where the optimum scheduling of corrections can be 
obtained analytically . 


The third assumption will be valid^ for the present error models for 
midcourse correction magnitudes to 5 m/sec. All but the first correction will 
usually fall within this range. Again^ this presents no difficulty since the 
timing of the first correction will be empirically determined. 


If the three assumptions are satisfied^ the optimum correction schedule 
is obtained as follows: Assimiing there are N corrections^ the optimum tim- 

ing of the jth correction is 


tF - tj . ^ 


niiy 


j = 2, 


N 


where mp 


is the allowable rms terminal miss and mx is the miss uncertainty 
after the first correction. The timing of the first correction is determined 
empirically such that the sum of the rms corrections 

N r ('+ 

, fm, ( t, ) 

T = ) Cj = Ci(ti) + (W - 1) TP I ^ ■ 


m. 


•a 


0 = 1 


IS a minimum. 


Lawden (ref. 2) reports an optimum number of corrections^ which give 

the least value of T. The value of N* is typically between 8 and 10. How- 
ever^ the present work will usually use a four-correction schedule since for a 
larger number of corrections the individual corrections become small^ of the 
order of the errors in making corrections^ and the time interval between the 
last few corrections also becomes small. Appendix A indicates that there is 
little fuel penalty for using four or five corrections rather than the optimum 
number . 
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RESULTS AM) DISCUSSION 


The following discussion considers the nature of reasonable observation 
schedules for the combined use of radar and on-board observations^ and com- 
pares the resulting performance with that theoretically predicted for the 
optimum correction schedule and with the perforaiance of the purely on-board 
system of reference 1. 


Observation Schedules 

The observation schedules are determined empirically. A schedule is 
sought that satisfies the condition that the miss can be better estimated than 
it can be corrected at the time of each correction. When this goal is 
approached, the fuel requirements will be near the minimum and the terminal 
miss performance can then be improved somewhat by increases in the number of 
observations. The initial choice of schedule and further adjustments to this 
schedule need not be made arbitrarily; past studies, such as references 1 
and 3 have already indicated that radar provides superior information in most 
regions of the missions and that on-board observations are most effective dur- 
ing approach and passage near the planets. Thus, an initial schedule might 
consist entirely of radar measurements. The number of measirrements is then 
increased, and on-board observations are added and increased where effective 
until further adjustments show very little Improvement in performance. 

More organized empirical approaches, such as the method of steepest 
descent employed in reference 4, can be used to generate the most economical 
schedule. However, economical schedules are not of interest here since ground 
tracking normally generates and processes data in large quantities. Rather, 
this study seeks sample schedules which demonstrate the performance possible 
when both on-board observational and radar -tracking data are used. 

Observation schedules for each leg of both reference missions are dis- 
cussed below. The observation and correction schedules are listed in table III, 
and the corresponding time histories of rms miss uncertainty are plotted in 
figures 3 and 4. It should be noted that, although these curves are drawn 
smoothly, the rms miss uncertainty is a decreasing step function with step 
changes occurring at each observation. In addition, accelerometer measurement 
errors introduce a cori 5 )aratively minute step increase in rms miss uncertainty 
at the time of each correction, although this effect is usually not visible at 
the scale' of the figures. 

High-sp eed mission. Earth -Mars leg . - The reference trajectory for this 
case is type”l~(trajectories will be classed as type I if the total heliocen- 
tric angle from injection to arrival is less than l 80 °, and type II otherwise) 
with a total flight time of 112.5 days. 

The navigation schedule (table Ill(a)) and the resulting time history of 
the miss uncertainty (fig- 3(a)) show that a modest number of radar observa- 
tions can rapidly reduce the miss uncertainty to several hundred kilometers. 
This accuracy will satisfy the minimum fuel condition until the second 


9 



INI III 


nil IH I ■■ 1 1 


II I IM ■ 


midcoiirse correction; therefore^ no observations are made during most of the 
heliocentric phase of the flight. During the approach to Mars^ radar observa- 
tions alone are inadequate; however, frequent theodolite observations give a 
sufficiently accurate determination of the miss. 

The range -rate measurements provide very little information on the miss 
due to the large measurement errors postulated. Simulation without range -rate 
measurements (results not shown) shows virtually the same performance as that 
obtained for the system with both range and range rate. Thus, the performance 
obtained here for the radar system would apply as well if only range measure- 
ments were made. Similar remarks apply to the remaining legs of the two 
missions . 

The approach -phase schedule can be examined in greater detail to illus- 
trate the capabilities of on-board and radar measurements in this region. The 
approach phase refers to the 3 -day interval from 109*3 to 112.32 days which 
includes the last two velocity corrections and ends at the final correction. 
For the reference navigation schedule, the rms miss uncertainty at the begin- 
ning of the approach phase is 36.6 km, which, after 42 theodolite observations 
of Mars, reduces to 4-9 km* 

The capabilities of theodolite observations during the approach phase are 
illustrated in figure 5* The data show the variation in the downrange, cross - 
range, and radial components of the terminal miss uncertainty for (l) differ- 
ent numbers of observations in the approach phase, and (2) three initial 
levels of miss uncertainty. The approach-phase schedules used here consist of 
a sequence of theodolite observations with twice as many between the third 
and fourth corrections as in the segment prior to the third correction. In 
each case the observation rate is uniform in each of the two segments of the 
approach phase. The minimum number of observations considered is 12 and the 
maximxmi 99^ "the corresponding maximum observation rate is once each half 

hour. The three levels of initial miss uncertainty resulted from the omission 
of various blocks of observations from the reference schedule prior to the 
approach phase. 

Two fimdamental differences are seen between the estimation of the down- 
range component of miss imcertainty and the crossrange and radial components. 
First, the performance in the downrange component is dependent upon and 
approximately proportional to the initial uncertainty; whereas, in the other 
components, performance is independent of the initial values. Second, the 
minimum number of observations gives nearly as good results for the downrange 
component as the maximum number, while in the radial and crossrange con^onents 
performance is significantly Improved when the number of observations is 
increased. This suggests an inherent limit on the ability of theodolite 
observations in this interval to resolve the downrange con 5 )onent of miss. 

To obtain comparable data on the usefulness of radar data during the 
approach phase, additional simulations were run (results not shown) which 
indicated that even very accurate radar observations in the approach phase 
(1 km in range and O .5 m/sec in range rate) cannot reduce the downrange, 
crossrange, and radial coi35)onents of terminal miss uncertainty much below the 
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initial values. Apparently^ for the particular reference trajectory used 
here^ none of these coirponents of miss uncertainty correspond very closely to 
that conponent which can be resolved by radar data. 

High-speed mission^ return leg . - The reference trajectory is type II with 
a flight time of 190.77 bays. 

With radar (table Ill(b)^ fig- 3(b)) no on-board observations are neces- 
sary. The schedule consists of radar measurements at 4-day intervals until 
the first correction and at 10 -day intervals thereafter until the last 15 days 
of the flight when more frequent observations are made. Relatively few obser- 
vations were required during the entire flight. 

It may be noted that observations of Mars during departure from that 
planet are omitted from the schedule. Reference 1 indicates that^ even with 
an entirely on-board system^ observations of Mars at this time are ineffective 
in estimating the miss compared to observations of the Earthy despite the 
large distance to the Earth. 

Since the reference trajectory is type II^ there is a guidance singular- 
ity^ which^ in this case^ occurs at 43 days. The first correction cannot be 
made near this time. For the present reference trajectory there is no per- 
formance advantage in making the first correction as early as possible after 
departing Mars. Thus^ the first correction occurs at 84 days^ providing the 
estimation system with a very long time to determine the miss. It should be 
noted that a somewhat different reference orbit may require an early first 
correction and^ hence^ an early deteimination of the miss. A case of this 
type occurs in the return leg of the Venus swing -by mission discussed below. 

Venus swing -by miss ion j Earth -Mars leg .- The reference trajectory is 
type I with a flight time of I 70 days. Table IIl(c) lists the navigation 
schedule used to generate the data for this leg^ and figure 4(a) shows the 
resulting miss uncertainty. The results are generally similar to those for 
the outbound leg of the high-speed mission and will not be considered further. 

Venus swing -by mission^ Mars-Venus leg .- The type I reference trajectory 
for this leg has a flight time of I 85 days and remains at distances of the 
order of 10® km from the Earth throughout the flight. The observation sched- 
ule (table 111 (d)) consists of some observations of Mars early in the flighty 
radar observations at l/2-day or 1-day intervals up to the time of the first 
correction and at 10 -day intervals thereafter^ and frequent observations of 
Venus during the final 20 days. Figure 4(b) shows the resulting time history 
of the rms miss -estimation error. 


4?he guidance singularity occurs because A 2 (t-p^t) becomes singular when 
the heliocentric angle to the target position is l80^. In this case^ the out- 
of-plane conponent of the miss is uncontrollable by velocity corrections and 
is poorly controllable in the region near the singularity. This does not pre- 
clude the possibility of correcting only the controllable conponents ^ but such 
sophistication of the guidance system is not within the scope of this study. 
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Radar measurements are very effective up to the approach to Venus ^ and 
are capable of reducing the miss -estimate error to several hundred kilometers. 
However^ prior to the first correction (at 30 days)^ the most effective sched- 
ule combined both Mars and radar observations. Figure J+(b) con 5 )ares the 
results for this region of the mission for three schedules: only on-board 

observational data^ only radar data^ and both radar and Mars observations . 

The schedules with radar data give a much better estimate of the miss than 
does the schedule with only on-board observations and the combined schedule 
does somewhat better than the radar -only schedule. This occurs because the 
radar is unable to estimate the crossrange miss during the early part of this 
leg. Figure 6 coTopSiTes the early time history of the error in estimating the 
crossrange miss for both the combined and radar -only schedules^ and shows the 
importance of the Mars observations for an early determination of this com- 
ponent. However^ in the present simulation^ this component is only a small 
fraction of the initial miss (cf. table II) ^ and the radar alone is able even- 
tually to make a good estimate of the crossrange miss. Thus^the principal 
effect of the early Mars observations is a small saving in the total velocity 
correction (about 1 m/sec). Of course^ if the initial crossrange miss were 
larger^ the Mars observations would give a correspondingly more significant 
fuel saving. 

Venus swing -by mission^ Ve nus -Earth leg.- The type II reference trajec- 
tory for this leg has a guidance singularity at 28 days^ and a total flight 
time of 125 days. The distributions of initial deviations and the imcertainty 
in estimating these deviations are obtained from simulation of the preceding 
Mars -Venus leg of the mission and are given in table II- 

The observation schedule (table Ill(e)^ fig- 4(c)) consists of frequent 
observations of Venus while in the immediate vicinity of that planet and radar 
observations throughout the flight. 

The fuel requirements for this leg are decreased if the first correction 
is made prior to the guidance singularity^ within the first 10 days. The most 
effective schedule for an early determination of the miss combines information 
from both the radar and observations of Venus. Figure 4(c) compares the time 
histories of miss uncertainty during the first 10 days of this leg for three 
schedules usings respectively^ only observations of Venus ^ only radar data^ 
and both radar data and observations of Venus. It is seen that the combined 
schedule does better than either of the schedules with only one type of obser- 
vation. Of course^ on the return leg the miss can be estimated very accu- 
rately with the radar alone and there would be no reduction in terminal miss 
performance if the Venus observations were deleted. The principal effect of 
the Venus observations in the combined schedule is then a modest reduction in 
velocity-correction requirements (about 3 m/sec in the present example). 

Conp)arison With Performance Predicted by the 
Optimum Correction -Schedule Theory 


The optimimi correction schedule and corresponding minimum fuel 
requirements are computed^ as described in appendix after an enipirical 
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determination of the optimum timing of the first correction. The timing of 
the final correction is determined by a specified rms allowable terminal miss 
which^ in the present case^ is 2 km. 

Performance data from simulation and from the theory are compared in 
tables IV and V for all legs of both reference missions. Data for three-^ 
four-^ and five -correction schedules for both legs of the high-speed mission 
are contained in table IV and data for four -correction schedules for all legs 
of the Venus swing-by mission in table V. The specified terminal miss per- 
formance is generally not achieved^ principally because the miss has not been 
determined with sufficient accuracy at the time of the final correction. How- 
ever^ performance is generally good because the largest component of the miss 
is oriented in the downrange direction. This affects only the arrival time as 
discussed later. 

In general^ the predicted fuel requirements are optimistic. This is 
because the theoretical model neglects the uncertainty in the miss estimate 
at the time of each correction and the miss generated by correction errors 
other than the cutoff error. Both of these approximations cause the size of 
the third and subsequent corrections to be underestimated by the theory. The 
total velocity correction shows good agreement for the four- and five- 
correction schedules^ being within 4 m/sec in every case. Table IV also 
illustrates the decrease in velocity-correction requirements and the size of 
each correction with an increase in the number of corrections. 

Thus,, because of the accuracy achieved with radar observations,, the 
theory of Lawden (ref. 2) gives an adequate first approximation to the optimum 
correction schedule^ velocity-correction requirements^ correction size^ and 
their dependence on the number of corrections. 

Comparison With Performance of the On-Board Navigation System 

The rms total velocity-correction requirements and terminal miss for the 
system of this report and a system using only on-board observations are com- 
pared in table VI. Data for the on-board system are taken from reference 1 
and correspond to schedules using four midcourse corrections on each leg with 
variable time -of -arrival (VTA) guidance corrections where indicated. Data for 
the present system correspond to the navigation schedules of table III and use 
fixed arrival -time guidance for all corrections. 

The results show a varying amount of fuel saving with the inclusion of 
ground-based radar data; for the return leg of the high-speed mission there is 
almost no reduction^ but for the outbound leg of the Venus swing -by mission 
the velocity-correction requirements are cut in half. 

Table VI shows a large reduction in downrange miss with the use of radar. 
Navigation data from on-board observations generally give limited information 
on the downrange miss until late in the mission,, where the fuel required to 
correct this miss is large. This deficiency of the on-board system made it 
advisable to utilize VTA guidance for the final one or two corrections in ref- 
erence 1 in order to reduce the fuel requirements substantially. This type of 
guidance corrects only the lateral miss and ignores the downrange miss 
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remaining from earlier fixed time -of -arrival corrections; hence it produces 
the large downrange miss seen in table VI for the on-board system. On the 
other hand^ the system of this study improves the estimate of all components 
of the miss everywhere in the flight; consequently^ only minor fuel saving 
•would occur if VTA corrections were made during the approach phases. 

The principal effect of the downrange miss is a change in the time of 
pericenter. This is seen from the following linearized equations for the 
errors in pericenter position and time (see ref. 5): 

= (ui? + nr?)6r(tp) - — pu^6v(tf) 

5tp = - ^ p^Sr(tf) - ^ 

where u^ p, and n form an orthonormal vector triad with u along the ref- 
erence pericenter position (radial) and p along the reference pericenter 
velocity (downrange). To first order, the downrange miss, ^6r(tf ) , affects 
only the time of pericenter. The reduction in downrange miss sho-wn in table VI 
corresponds to a reduction in the rms pericenter time error from several min- 
utes for the on-board system to several seconds. However, the timing error 
associated with the on-board system requires only modest adjustments in entry 
range for landing at a fixed site on Earth. 

Of particular interest in manned missions is the control of pericenter 
radius to lie within a safe entry corridor both at Mars and at the critical 
retixm to Earth. It can be seen from the above equation that the error in 
pericenter radius, u^&Fp, is simply the radial miss, which, in table VI, is 
generally cut in half when radar data are Included. 


CONCLUDING REMARKS 


This paper examines the performance of midcourse guidance and navigation 
systems for manned interplanetary missions when both radar and on-board obser- 
vational data are available. 

The investigation of appropriate observation schedules for such systems 
shows that radar generally provides superior data and will be sufficient for 
best performance during most of the mission. However, during the approach to 
distant planets on-board observations are essential to Improve the miss esti- 
mate to the accuracy required for entry or for close approach of manned 
vehicles. In addition, on-board observations during departure from distant 
planets will usually reduce fuel requirements. 

It may be noted that the amount of radar data assumed in this study is 
quite small compared to what would be available in the real case. More inten- 
sive use of radar data would improve performance somewhat over that shown in 




this report, but the iirprovement would be limited because the trajectory is 
already estimated with the minimal data nearly as accurately as it can be 
controlled. 

When the trajectory is estimated more accurately than it can be 
controlled, as in the present study, the optimum correction schedule is inde- 
pendent of the observation schedule. This feature, together with several sim- 
plifying assumptions, allows application of the work of Lawden (ref. 2 ) to the 
a priori coirputation of an optimimi correction schedule. The minimum velocity 
correction requirements, correction size, and the dependence of these param- 
eters on the number of corrections are readily determined by the theory. An 
optimum number of corrections, on the order of 8 to 10 for the cases studied 
here, is given by the theory. However, this can be reduced to 4 or 5 cor- 
rections without significant fuel penalty. Results from simulation of the 
system gave fair agreement with the theoiy. 

The combined radar and on-board data are shown to give fair improvement 
in both fuel requirements and miss performance over a purely on-board system. 
However, it should be noted that the on-board system for which the comparison 
was made used a minimal amount of data and sufficient additional data could 
be expected to improve its performance. 

The performance obtained from simulation of the guidance and navigation 
system of this study gave velocity-correction requirements ranging from l8 to 
24 m/sec on the various legs of the two missions. The rms miss in pericenter 
radius was less than 2.8 km at Earth and Mars and 6 km at Venus. The rms 
downrange miss varied from 2 to 22 km. 


Ames Research Center 

National Aeronautics and Space Administration 
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APPENDIX A 


DERIVATION OF OPTIMUM VELOCITY- CORRECTION SCHEDULE 


The correction schedule is designed to null the miss. The optimum 
schedule is one which accomplishes this end at least cost^ that is^ for the 
least amount of corrective fuel. In general^ the fuel requirements vary with 
the observation schedule as well as the correction schedule^ and the general 
problem is of such complexity that no analytical solution has thus far been 
given. However^ because of the inclusion of radar data in this study ^ the 
miss is generally estimated better than it can be controlled. Consequently^ 
certain simplifications can be made. Firsts the optimum correction schedule 
will^ for the most part^ be independent of the observation schedule. Secondly^ 
the first midcourse correction will null almost all the miss generated by 
injection errors^ and subsequent midcourse corrections will correct principally 
for miss generated by random errors in the preceding correction. Of course^ if 
the estimation and correction systems are so accurate that the miss following 
the first correction is acceptable^ then the scheduling problem will be tri- 
vial. But for the assumed state of the art^ the miss generated by errors in 
each correction is significant. This miss has a lower bound determined by the 
thrust cutoff error and is approximately proportional to the time to go. 

In the following^ the general equations for the statistical properties of 
the velocity corrections are derived and a statement of the optimum scheduling 
problem is given. Next^ a list of simplifying assumptions is discussed and^ 
with these assuir^tions^ an optimum correction schedule is derived. The sim- 
plified problem has been solved by Lawden (ref. 2)^ and the principal features 
of that solution and its application to the present simulations will be 
discussed. 


General Problem 

Velocity corrections.- The guidance correction is an attempt to null the 
miss. For fixed time -of -arrival (FTA) guidance^ the miss is the position 
deviation from the reference position at the reference time of arrival^ tp.. 

If the quantities of interest are assumed to be adequately represented by the 
linearized equations of motion^ the miss is related linearly to the current 
state deviation: 

m(t) = B(tp^t)x(t) (A1) 

where ^ for FTA guidance^ B is the 3x6 matrix taken from the upper half of 
the transition matrix (ref. 1): 

B(tF^t) = [Ai(tp^t) A2(tpi^t)] (A2) 

More generally^ for arbitrary linear guidance laws^ the miss can be expressed 
as in equation (Al)^ but with different expressions for B (ref. 5). 

The current state^ x^ can be given in terms of its estimated value^ x^ 
and the error in the estimation^ x^ as 

l6 



similarly^ 


(A3) 


where 


x(t) = x(t) + 5c(t) 

m(t) = m(t) + m(t) 

m(t) = B(tp,t)x(t) 
m(t) = B(tp,t)x(t) ^ 




The conmianded velocity correction is con 5 )uted from the estimated miss . 
(a justification for this procedure is given in reference 6.) For FTA guid- 
ance^ the correction is 


AVj, = -A2^(tp,t)m(f ) (A5) 

where the symbol t" refers to time ^just prior to t. This velocity correc- 
tion is carried out with random error, 3 ^ so that the actual correction is 

AV = AVc + 3 (a6) 

The miss vector . - The actual miss is a step function with changes at each 
correction. Its initial value is specified by the injection errors, XqJ 
hence, prior to the first correction, the miss is 

m(t) = B(tp,tQ)xQ, to < t < tj (A7) 

After each correction, the miss is given by the miss prior to the correction 
and the change due to the correction. In the interval (tj,tj^^) between the 
jth and J+lst correction, the miss is: 

m(t) = m(tj) + A 2 (tp,tj)AVj 

From equations (A3) and (a 6), this becomes 

m(t) = m(tp + A 2 (tj,,tj)rij , t^ < t < t^+i (A8) 

That is, the miss following a correction is due to the miss uncertainty prior 
to the correction and the error in making the correction. 

S tatistical properties . - The navigation and guidance process is subject 
to random errors^ so that the parameters which give the performance of the 
system - the final-miss and total -fuel - are random variables that must be 
dealt with statistically. All input random variables - injection^ measure- 
ment^ and correction errors - are assumed to be gaussian distributed with 
zero means. Since the equations of motion have been linearized^ the random 
vectors of interest are linear combinations of the input random variables and 
are^ therefore^ gaussian distributed with zero means. In that case^ their 
statistical properties are specified by their covariance matrices. 
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AssiJiming the state estimate has "been made using the Kalman filter^ then 
X and X are uncorrelated and^ consequently^ m and m are also uncorrelated. 
The covariances and M are defined by: 

M = E[imiiT] y etc . 

where E is the expectation operator. These are positive definite matrices 
related by 


M = M + M 


(A9) 


From equations (A7) and (a 8)^ one can obtain 

M(t) = B(tF,t)PARoB^(tp,t) = Mq , 
M(t.) = + A2('tp';'t)SjA2(tjp^'b) f 


tQ < t < ti 

< t < + 


(AlO) 


where PARq Sj are the covariance matrices, respectively, of the injec- 
tion errors and the errors in making the jth correction. The second of 
these equations uses the independence between correction errors, -qj, and 
estimation errors, m(tj). 


Define the covariance matrix, Cj, as E[AVq AVq] . This can be written, 
from equations (A5) and (A9), as 


Cj = A|i(tF,tj)[M(tJ) -M(tJ)]A2^ (tp,tj) 


(All) 


•The actual correction (eq. (a 6)) then has the covariance 


Vj = ^j ^ 


(A12) 


where statistical independence between the commanded correction and the cor- 
rection error has been assumed. 

Problem of the optimum correction schedule . - If N midcourse velocity 
corrections are made between injection^ to^ and arrival^ tp^ their execution 
times^ t-^] y are called the correction schedule. The optimum 

scheduling problem to be discussed is stated as follows : 

Given that N corrections are made in (tQ^tp)^find the sequence of 
execution times 


t^ < ti < . . • ^ ^ *^F 

such that the sum of the rms velocity corrections ^ 


T = 


N 

^/trace(Vj) 

J=i 
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is a minimum,, and the rms miss following the last correction is 




= 


Here^ is some acceptable n-umerical value of the rms final miss. 

The parameter mg, represents the desired terminal miss performance and 
affects velocity -correction requirements. A lower value of mg^ (i.e.,, better 
terminal miss performance) increases the required total correction; that is^ 
the choice of m^ involves a trade-off between fuel and terminal miss. It 
should also be recognized that^ regardless of the choice of m^,, the terminal- 
miss performance can be limited by other factors inherent in the navigation 
system. 

The choice of cost function T^ as given above,, is motivated by the need 
for mathematical tractability . A more appropriate criterion would be the 
minimization of expected fuel requirements which,, for a vehicle of constant 
mass^ is given by the sum of the expected velocity -correction magnitudes 
rather than the sum of rms values : 

T* =^E {jAVjl} 

J=1 

However, if the magnitudes of the several corrections are similarly distri- 
buted (i.e., the distribution curves have the same shape and differ only in 
scale), then the rms and mean magnitudes of each correction are always in the 
same ratio. Previous experience shows that such proportionality is a fair 
approximation in actual practice. Thus, minimizing the sum of rms correction 
magnitudes is sufficiently close to minimizing the fuel requirements to Jus- 
tify the choice of the former as the optimality criterion. 


Simplifying Assumptions 

Independence of correction and observation schedules .- For convenience, 
equations (AlO), (All), and (A12) are repeated: 

For J = 1, 2, . . . , N 

Vj = A2^(tp,tj)[M(tj) - M(tj)]A2^^(tp,tj) + Sj 


M(t) 


Mq tQ < t < tl 

M(tj) + A2(tp,tj)SjAl(tj.,tj) tj < t < 


In general, the total fuel depends on the performance of the estimation system 
through the terms M(tj), J = 1, 2, . . . , N. The analytical problem pre- 
sented by the general case is unsolved, but adequate schedules can be obtained 
enpirically (e.g., ref. U) . 
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On the other hand^ if the term M(tj) can be neglected^ then the correc- 
tion schedule is independent of the observation schedule for the interval 
(tj.j^^tj). Further, if all the terms M(tj) j = 1, . . . , N can be neglected, 
the correction schedule is entirely independent of the observation schedule. 

It is clear that if the miss were perfectly estimated at the time of the 
Jth correction then M(tj) would be zero. In general, M(tJ) is not zero and 
it is of interest to state how well the miss must be estimated in order to 
neglect M(tJ). 

Assumption 1 . 

stt" T) can be neglected if 
J 

(a) tr [M(tJ)] < < tr [M(tj)] 

(b) tr [M(tJ)] < < tr [A 2 (tp.,t j )SjA 2 (tp,t j ) ] 


Ibe first condition states that the miss uncertainty is small conpared to the 
miss, from which 

Vj = A2^(tp,tj)M(tT)[A2'^(tj.,t^)]^ + Sj (A13) 

The second condition states that the miss has been estimated better than it 
can be corrected, that is, the miss following the jth correction is due 
almost entirely to errors in making the correction. 


M(t) 


Mq t^ < t < ti 

A 2 (tj.,tj)SjAl(tp,tj) tj < t < t^+^ 


(Al4) 


Simplified correction -err or model.- The covariance matrix, Sj, is an 
important parameter when the correction and observation schedules are inde- 
pendent. The error model, already discussed in the text, is taken from ref- 
erence 1, appendix A. 








C. 


where 7 ^, and are the variances, respectively, of the proportional, 

aiming, and cutoff errors, and c^ is the mean-square correction, tr(Cj). 




If the commanded correction is almost spherically distributed so that 
1 ^ T _ +.Vipn s can he approximated by 


T c 2 I , then 

3 3 


"a " 




+ 7^)c5 


+ e" 


ii 
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For the present Mars mission studies, the rms errors are assumed to be 



0 . 01 , 





= 0.2 m/ 


sec 


i — In this case^ the cutoff error is the dominant error ^ statistically^ when 
>7cj is about 5 m/sec or less. The first correction is usually large^ statis- 
tically^ but the remaining corrections are usually small, allowing the follow- 
ing approximation: 


Assijmption 2 . 

For each correction except the first, the commanded velocity 
correction is distributed almost spherically and the cutoff error 
is, statistically, the dominant error in making the correction: 



j = 2, 


N 


(A15) 


Approximation of Vj— The actual correction consists of the commanded 
velocity correction and the errors in executing the correction (eq. (A12)). 
The error term is usually small by comparison, so that one can assume: 


Assumption 3 > 

The errors in making a correction are small co 2 i 5 >ared to the 
actual correction 

tr(Sj) « 

In this case, equation (A12) becomes 




(Al6) 


The first correction is usually large and assumption 3 readily applies. 
For the remaining corrections, in view of assumption 2, assumption 3 is satis- 
fied if « c2. For the assu med value nF = 0.2 m/sec, assumption 3 is 
valid when the rms correction, is 1 m/sec or greater. This should nor- 


mally be the case, since the operation of thrust devices near the level of the 
errors in thrusting is generally avoided. 


Approximation of A^Ctp,!).- The columns of this matrix are the deriva- 
tives of position at tp with respect to the components of velocity at t, 
evaluated on the reference orbit. 


If the gravitational acceleration on the deviated path is identical to 
the gravitational acceleration on the reference path at every time, then the 
deviation equations of motion are those of a straight line; therefore 


^sC'bp^t) (tp - t)I 


(A17) 


For the heliocentric orbits occurring in interplanetary missions, the differ- 
ence in gravitational acceleration on nearby orbits is very small so that 
equation (A17 ) remains a fair approximation over extended periods of time. 
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The maximunL interval for which this approximation remains valid is not pre- 
cisely known^ but an absolute upper limit is given by the fact that AsCtp^^t) 
becomes singular in Keplerian motion when the interval tp-t corresponds to 
a difference of l 80 *^ in true anomaly. 

Equation (AI7) is used in the remainder of this analysis in the equations 
governing the second and later corrections^ which normally occur when much 
less than half an orbit remains to be traversed to the target position. How- 
ever^ equation (AI7) is often a poor approximation at the time of the first 
correction; The matter of approximations to the transition matrix is dis- 
cussed at length in reference 7 • 


Simplified Problem and Its Solution 

Qptimim correction s chedul e for the simplified problem . - In the simpli- 
fied problem only assumption 3 is applied for the first correction. All three 
assumptions and the approximation given in equation (AI7) are used for the 
remaining corrections . 

If ci(ti)^ nii(ti) are the rms first correction and miss following the 
first correction^ then 


Vtr(Vi) = ci(ti) 

ArCVs) = mi(ti)/(tF-t2) 

N/tr(Vj) = 


j = 3 , 


,, N 


The dimensionless time to go is defined as 


Tj = (tp - tj)/(tF - to) . 


The sum of the rms correctioas is then: 


T = ci(ti) + 




1 . ^ 


tp-to 


N 


(AI8) 


J =3 


From equations (Al 4 ), (A 15 ), and (AI7) the rms final miss is entirely deter- 
mined by the cutoff error in the last midcourse correction, 


m(t^) = (tp - 


Since the rms final miss is specified to be equal to 1%, then the timing of 
the final correction is specified by 


Tm = 


“a 




(AI9) 
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It should be noted that the sin^lified equation for the final miss accounts 
only for the miss generated by the cutoff error in the final correction and 
therefore underestimates the true value of the final miss to the extent that 
the miss estimate uncertainty at tjj contributes to the terminal miss . 


The independent parameters of the correction schedule are now ti^ 
^ extreme (minimum) at 


5T 

Sti 


= 0 


= 0 


From equation (A18) 




9ti Sti tTp-t2 5ti 


j = 2, 


= 0 


N-1 




ST 

dx2 




■fep-to 




/ (A20) 


— = n/? (- =0 0 = 3, . . N-1 

St. V t 2 


j+1' 

“ u 

Except for the firsts these equations are readily solved (algebra omitted)^ 
giving j -1 


N-i 


m. 






j = 2, 


N 


(A21) 


That is^ the optimum schedule for ts^ • • • ^ depends on the timing of the 

first correction only through the rms miss^ mi^ inherited from the first 
correction. The sum of the rms velocity corrections is now 

1 


T = cx + 7^(N-1)(^^ 


N-1 


(A22) 


where T is given as a function of ti only. If ti can be determined 
(empirically) such that T is a minimum^ then all of equations (A20) are 
satisfied and the entire schedule [ti^ . . t-^] is the solution of the 

optimum scheduling problem. 


Several consequences of equation 
the jth correction is a fixed fraction 
correction: ^ 



should be noted. The timing of 
the timing of the preceding 


2. 


N 


(A23) 
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For all corrections after the firsts the rms miss is reduced by a fixed 
fraction: 



and the rms correction is constant: 


1 



(A24) 


(A25) 


Equation (A18) and its solution are given in reference 2. The optimum 
schedule represents a compromise between two opposing effects. If the jth 
correction is delayed^ the miss generated by random errors in the jth cor- 
rection will be less and subsequent corrections will require less fuel. On 
the other hand^ a delay increases the fuel required for the jth correction. 


Computation of tx > - The timing of the first corrections is left to 
empirical determination. Hence^ only assumption 3 need apply at tj,. 

It is not difficult to determine ti empirically if one assumes a 
reasonable observation schedule for the initial part of the flight and then 
uses a conputer program which simulates the statistical performance of the 
guidance and navigation system^ obtaining mi(t'*") and Ci(t). These parameters 
are^ respectively^ the rms postcorrection miss and velocity correction if a 
correction is made at t. Equation (A22) can be used to conpute the total 
velocity-correction requirements as a function of the time of the first 
correction^ from which the optimum time^ ti^ is selected. 

Number of corrections . - An optimum value of N^ found by minimizing T 
(eq. ( 22) ) with resp^^t ^fco is reported in reference 2. The optimum N 
is one of the neighboring integers of 


N* = 1 + Z 

For typical values of the ratio mi/mQ, found in this study ^ 10^ to 10*^, we 
have 



N* = 8.1 to 10.2 

However^ practical difficulties^ which will be discussed below^ arise with 
such a large number of corrections^ and the fuel saved by using 8 or 10 cor- 
rections rather than^ say^ 4 or 5 may not be significant. 

Data illustrating the dependence of fuel requirements and correction 
magnitudes on the number of corrections are given in sketches (a) and (b). 

The first correction^ which accounts principally for injection errors^ is very 
nearly independent of the number of corrections^ but fuel requirements for the 
remaining corrections depend on N (see eq. (A22)). Sketch (a) gives the sum 
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Sketch (a) 



Sketch (h) 


of the rms velocity corrections for all 
corrections except the first versus the 
parameter mi/ma for various values of 
N. It is seen that the fuel required 
decreases with increases in the number 
of corrections^ but for large N the 
advantage of further increases in N 
diminishes . 

Sketch (b) shows that the size of 
individual corrections also decreases 
with increasing N and soon approaches 
the size of the cutoff error. Note 
also that for the low values of N the 
individual corrections can be large 
enough that significant miss may be 
generated by correction errors other 
than the cutoff error. In this case^ 
poorer agreement with the theoretical 
fuel requirements is expected. 

A final point of interest is that 
the time interval between the last few 
corrections becomes small as N 
increases. This can be seen from equa- 
tion (A23 ) . 

In the simulations of this study ^ 

N = 4 was selected as a compromise 
between total fuel and correction size. 

Off -optimum timing . - It is of 
interest to determine how critical the 
timing is of each correction. 

Usually, the fuel requirement per 
unit miss rises very slowly with time 
in the region of the first correction 
(although this is not always the case). 
In the present simulation, the miss is 
determined at the time of the first 
correction to the order of the miss 
generated by errors in the first cor- 
rection or better. In this case, there 
is very little penalty associated with 
delays in the first correction. In 
some cases, a guidance singularity 
occurs soon after injection so that 
delays beyond the optimum timing of the 
first correction may cause more sig- 
nificant penalties. 
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On the other hand^ if the first correction is made earlier than the 
optimum time, the principal penalty arises from a poorer estimation of the 
miss. For instance, if the correction is made early enough that the post- 
correction miss is twice the value for the optimum timing, the fuel for the 

remaining corrections will be increased by a factor of while the 

fuel for the first corrections will remain very nearly the same; that is, 

T(2mi) = T(mi) + - 1] [T(mi) - Ci] 

where mi refers to the miss following the optimally timed first correction. 

For corrections 2, . • N-1 (t]\j- remains fixed to satisfy the miss 

performance criterion) , let the timing of the jth correction be changed from 
the optimum time, tj, by the amount At^ . This yields the change in the sum 
of rms velocity corrections (algebra omitted) 

AT _ [Atl/(tF-tj)] 

In this equation the penalty, AP, has been made dimensionless by a factor 
recognized as the rms size of the individual corrections. The penalty is 
greater for a delay (positive Atj) than for correcting early by the same 
amount, and, of course, it becomes axbitrarily large as Atj approaches 
tp-tj. A penalty of 1 percent of the individual correction is obtained for an 
off -optimum timing of 


tp-tj 


± 0.1 


Similarly, a 10-percent penalty is incurred for timing errors of 

At^ 

- +0.27, -0.37 

J 


This is considered a small penalty for fairly large percentage changes in 
timing and the correction schedule appears to have a broad optimum. 
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TABLE I.- REFERENCE TRAJECTORY PAilAMETERS 



Departure conditions 

Arrival conditions 

Trajectory 
and leg 

Date 

Velocity^ 

km/sec 

Altitude^ 

km 

Trip time, 
days 

Date 

Velocity, 

km/sec 

Altitude, 

km 

High speed 
Earth to Mars 

May 31.57, 1971 

11.761 

160.00 

112.43 

Sept. 21 . 00 , 1971 

8.49 

-3.54 

High speed 
Jfers to Earth 

Sept. 27.22, 1971 

9.813 

300.32 

190.77 

April 4.98, 1972 

14.26 

-2.96 

Venus swing -by- 
Earth to Mars 

Sept. 9.49, 1975 

12.020 

159.27 

170.01 

Feb. 26.49, 1976 

7.76 

26.43 

Venus swing -by 
Mars to Venus 

]Vte.rch 27.55, 1976 

7.097 

499.90 

185.54 

Sept. 29.10, 1976 

14.69 

3363.59 

Venus swing -by 
Venus to Earth 

Sept. 28.54, 1976 

14.992 

3349.96 

124.94 

Jan. 31 . 48 , 1977 

13.91 

-11.06 


Velocity and altitude are at periapse of the vacuum hyperbola. Dates are given in universal time. 








TABLE II.- RMS INJECTION^ CORRECTION, AND MEASUREMENT ERRORS^ 


Injection errors (Mars and Earth injection)^ 
Altitude 3.2187 km (2 miles) 4.47 m/sec (lO mph) 

Downrange 4.8285 km (3 miles) I.788 m/sec (4 mph) 

Crossrange 1.60935 km (l mile) 1.341 m/sec (3 n^h) 


Corresponding initial miss, km 


Mission Leg 

Altitude 

Downrange 

Cross range 

Total 

High speed Outbound 

62^000 

170,000 

27,600 

183,000 

Return 

104^000 

114,000 

9,200 

155,000 

Venus swing “by Outbound 

193,000 

198,000 

79,400 

288,000 

Mars -Venus 

122,000 

107,000 

5,940 

163,000 

Initial miss^- 

return leg. 

Venus swing -by 

mission, km 



Altitude 

Downrange 

Cross range 

Total 

Miss 

187,000 

172,000 

8,260 

255,000 

Miss uncertainty 

164,000 

146,000 

2,840 

219,000 


Errors _in making velocity corrections 

Magnitude, n/^ 1 percent 

Direction, 1^ 

Cutoff, 0.2 m/sec 

Errors in measuring velocity corrections 
1 cm/sec rms equally likely in all directions 
Theodolite observation errors 

10 seconds of arc 
0.001 
1 km 


Instrument error, \/^j_nst 
Radius uncertainty /pla net r adius. 
Position uncertainty^ v5R; 


4 


^All errors are assumed gaussian distributed with zero means. 

2 

The initial uncertainty is identical to the injection errors except on 
the return leg of the Venus swing -by mission. 

^These values are obtained from simulation of the preceding Mars -Venus 
leg, using the navigation schedule of table Ill(d). 
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TABLE III.- NAVIGATION SCHEDULE^ 
(a) High-speed mission, oufbo\md leg 


Initial Observation 

t ime y in t e r val y 

days days 


Group 

interval^ 

days 


Observation 

Number Group 


Type 

(2) 


0.1 

0.2 

5 

1.10 


Velocity correction 

1.50 

.5 

2.0 15 

93-50 


Velocity correction 

94.00 

.1 

1.0 75 

109.0 

.1 

2 

109.3 

.123 

14 

110.95 


Velocity correction 

111.00 

.048 

28 

112.32 


Velocity correction 

112.432 


Pericenter (Mars) 


Total observations- 42 radar^ 10 Earthy 87 Mars 


R 


R^2E 


2R,3M 

R 

M 


M 


^The table is read as follows: line 3 ^ for example^ specifies that each 

group of observations consists of one radar and two theodolite observations 
separated by 0 . 5 -day intervals. This group is repeated 10 times, beginning 
every second day. 

^R - radar; E, M, V - theodolite observations of Earth, Mars, or Venus. 
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TABLE III.- NAVIGATION SCHEDULE^- Continued 


(Td) High-speed mission, retvirn leg 


Initial 

Observation 

Observation 


time, 

days 

interval, 

days 

Nixniber 

(2) 

3-0 

4.0 

20 

R 

81).. 0 


Velocity correction 


101.0 

10.0 

8 

R 

172.0 


Velocity correction 


174.0 

2.0 

8 

R 

189.0 

.1 

3 

R 

189.3 


Velocity correction 


189.4 

.2 

6 

R 

189.5 

.05 

3 

R 

190.65 


Velocity correction 


190.765 


Perigee 



Total - 48 radar observations 



footnotes page 30. 
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TABLE III.- NAVIGATION SCHEDULE^- Continued 
(c) Venus swing -by mission^ Earth -Mars leg 


Initial 

time^ 

days 

Observation 

interval^ 

days 

Group Observation 

interval^ 

days Number Group 

.(2). 

0.1 

0.1 

8 

R 

•9 

Velocity correction 


1.0 

•5 

5 

R 

3.0 

1.0 

7 

R 

100.0 

10.0 

5 

R 

1U3.6 

Velocity correction 


145.0 

1.0 

4 16 4 

R, 3 M 

165.0 

1.0 

7 

M 

168.76 

Velocity correction 


168.28 

.02 

31 

M 

169.0 

.02 

44 

M 

169.90 

Velocity correction 


170.008 

Pericenter (Mars) 



Total - 29 radar ^ 9k Mars 


^•^^See footnotes page 30. 
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TABLE III.- NAVIGATION SCHEDULE^- Continued 
(d) Venus swing -by. Mars -Venus leg 


Initial 

Observation 

Group 


Observation 


time. 

interval^ 

interval^ 




days 

days 

days Nimber 

Group 

T^pe 

(2) 

1.0 

0.5 


4 


R 

S -0 

.5 

2 

30 

10 

2 R,M 

22.5 



1 


R 

23.0 

1 

2 

4 

2 

M,R 

27.0 


Velocity correction 




30.0 

1 

10 

l 4 

7 

R,M 

100.0 

1 

10 

12 

6 

R,V 

160.0 





R 

160.1 


Velocity correction 




160.5 

• 1 


5 


V 

161.0 

.5 


16 


V 

170.0 

.2 


1 


R 




63 


V 




1 


R 

183.85 


Velocity correction 




183.9 

.05 


SO 


V 

185.4 


Velocity correction 




185 . 45 



1 


V 

185.545 


Pericenter (Venus) 





Total - 43 radar, I 9 Mars, 121 Venus 
lessee footnotes page BO- 


SS 



TABLE III.- NAVIGATION SCHEDULE^- Concluded 
(e) Venus swing -by mission 


Initial 

time^ 

days 

Observation 

interval, 

days 

Group 

interval^ 

days Number 

Observation 

Group 


0.2 

0.02 

30 


V 

.62 

.03 

2 


R 

• 7 

• 1 

13 


3V,R,9V 

2.0 


Velocity correction 



2.5 


1 


R 

3.0 

1.0 

2.0 10 

5 

V,R 

20.0 

10.0 

8 


R 

95.0 


1 


R 

97.75 


Velocity correction 



100.5 

1.5 

15 


R 

123.0 


1 


R 

123.17 


Velocity correction 



123.2 

.2 

9 


R 

124.83 


Velocity correction 



124.945 


Perigee 




Total - k9 Venus ^ 42 radar 




^^^See footnotes page 30 * 
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TABLE IV.- COMPARISON WITH PREDICTED PERFORMANCE (RMS) 


(a) High-speed mission, outbound leg; flight time = 112.432 days 


N = 3 


N = 4 


N = 5 


Theory Simulation 


Theory Simulation 


0.55 

107.00 

112.32 


Total AV, 
m/sec 

Terminal 
miss, km 


10.77 

9-20 

9.20 


10.77 

8.93 

17.00 


29.17 36.72 


1.10 11.03 11.03 

93.50 2.56 2.42 

110.95 2.56 3-60 

112.32 2.56 3-41 


18.70 20.48 

2.0 5-4 


Theory Simulation 

0.24 10.37 10.37 

71.00 1.38 1.19 

106.50 1.38 1.46 

111.60 1.38 1.45 

112.32 1.38 1.66 

15.90 16.19 

2.0 3-4 


(b) High-speed mission, return leg; flight time = 


2.0 

190.765 days 


84.0 

185.5 

190.65 

10.15 

9.10 

9.10 

10.15 

7.94 

12.43 

84.0 

172.0 

189.30 

190.65 

10.15 

2.55 

2.55 

2.55 

10.15 

2.28 

2.67 

5.36 

84.0 

155.30 

185.50 

190.00 

190.65 

10.15 

1.35 

1.35 

1.35 

1.35 

10.15 

1.19 

1.41 

1.86 

1.64 

Total AV, 
m/sec 

28.35 

30.52 


17.80 

20.48 


15.55 

16.30 

Terminal 
miss^ km 

2.0 

4.5 


2.0 

3.8 

1 

2.0 

3.1 


Note : N = number of corrections 

ti = time of Jth correction 

c. = rms value of jth correction 
J 




































TABLE VI.- COMPARISON OF PERFORMANCE WITH ON-BOARD SYSTEM^ 


Mission 

High speed 

1 Venus swing -by 1 

Leg ; 

Outbound 

Return 

; Outbound 

Mars - Venus Return 

Navigation 

data 

On- 

board 

(2) 

On-board 
and radar 

On- 

board 

(2) 

Radar 

On- 

board 

(2) 

On-board 
and radar 

On- 

board 

(3) 

On-board 
and radar 

On- 

board 

(3) 

On-board 
and radar 

Total AV 

28.86 

20.48 

20.84 

^ 20.84 

41.61 

20.82 

i 3^-30 

18.37 

29.75 

23.45 

Radial ; 

miss ; 

4.8 

2.8 

5.0 

2.2 

4.6 

2.1 

9-6 

5.6 

4.2 

1.8 

Downrange i 
miss 

1173 

4.0 

741 

2.2 

590 

j 

17.4 

540 

21.3 

518 

2.1 

Total 

miss 


5.3 


3.8 


17.6 


22.1 


3.2 


^Four midcourse corrections on each leg 

^Last tvo corrections use variable arrival -time guidance 
^Last correction uses variable arrival -time guidance 








RMS miss uncertainty, 10^ km 


Observations 



Flight time, days 

(a) Outbound leg 

Figure 3*- MS miss estimate error - high-speed mission. 


H 





RMS miss uncertainty, 10 km 


Observations 

j^RRRRflR R R R R R R 
8R 


Initial value 



Flight time, c 

(a) Outbound le 
Figure i|-.- RMS miss estimate error, 


R 


R RyR Rj^R h ^ Hi* —j 
3M 3M lOM 75M 


M Mars observation 
R Radar measurement 
i Velocity correction 



Venus swing -by mission. 


Observations 


H RV RV RV RV RV RV RV RV RV RV RV 

5R MRR. MRR, - MRR 



(b) Mars -Venus leg 
Figure 4 .- Continued. 


RV 


RV 


Rl 


R 


II4V 


. Mars observation 
» Venus observation 
« Radar measurement 
I Velocity correction 





RMS miss uncertainty, 10 km 












Miss uncertainty at start of approach phase: 



Total 

Radial 

Crossrange 

Downrange 

o 

36.6 

21.0 

16.4 

25.1 

A 

100.4 

50.4 

47.0 

73.1 

□ 

158.2 

75.8 

77.6 

1 15.2 


-Downrange miss 


S --=8 


Radial miss 
Crossrange miss 

I 


Number of theodolite observations 


Miss performance variation with approach phase 
schedule, high-speed mission, outbound leg. 



RMS crossrange miss uncertainty, kn 


Schedules 

Radar-only 


R R 

RR 

RR 

R R 

R R 

R R 

R R 

R R 

RR R 

R RM 

R RM 

RRM 

RRM 

RRM 

RRM 

RRM 

RRM 




Radar only schedule 


•Reference schedule 


5 10 15 20 25 

Flight time, days 

Figiire 6.- Mars -Venus leg, effect of Mars observations on 
crossrange miss uncertainty- 


NASA- Langley, 1966 
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^^The aeronautical ayid space activities of the United States shall be 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phe 72 omena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof 

— National Aeronautics and Space Act of 1958 
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